/* COMPUTATIONAL FLUID DYNAMICS FOR AEROSPACE 
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ASSIGNMENT - I 
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USING C LANGUAGE*/ 
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c- PROGRAM 


/* AREA MACH NUMBER RELATION */ 


#include<stdio . h> 
#include<math . h> 


int main ( ) 

{ 

int i ; 

float M, G, A, B, C, D, E; 
for (i=l; i<=6; i++) 

{ 

printf ( "Enter the value of Mach number and specific heat ratio ") 
scanf ("%f%f", &M, &G) ; 

A=l+ (G-l) *M*M/2 ; 

B= (2/ (G+l ) ) *A; 

C= (G+l) / (2* (G-l) ) ; 

D=pow (B, C) ; 

E=D/M; 

printf ("The Area ratio A/A*=%f \n", E) ; 

} 

return 0; 

} 

OUTPUT: 


"F:\C works\Area mach.exe" 


ource» UIUEKKUMAR P» 

AEROSPACE TECH» MIT, MADRAS 

he Area ratio A/A* for 
1.05 1.40 

various Mach number <M> and specific heat ratio<G> 

1-0 . 050000 
1.10 1.40 

G =1.400000 

A/A*=ll. 591441 

1=0.100000 
1.15 1.40 

G =1.400000 

A/A*=5. 821828 

1=0.150000 
1.05 1.33 

G =1.400000 

A/A*=3. 910343 

1=0.050000 
1.10 1.33 

G=l. 330000 

A/A*=ll. 681849 

1=0.100000 
1.15 1.33 

G=l. 330000 

A/A*=5. 866473 

1=0 .150000 

G=l. 330000 

A/A *=3 .939487 

*ress any key to continue . . . 



/* PROGRAM TO FIND SUBSONIC & SUPER SONIC MACH NUMBER AT A 

PARTICULAR AREA RATIO*/ 


#include<stdio . h> 

#include<conio . h> 

#include<math . h> 

int main ( ) 

{ 

int i,j; 

printf ( " * AREA - MACH NUMBER RELATION *\n\n"); 

printf ( "Decide the flow condition \n\n"); 
for ( j=l; j<=6; j++) 


printf ("1. To find Area ratio at subsonic region\n"); 
printf ("2. To find Area ratio at supersonic region\n"); 
scanf ( "%d\n\n" , &i ) ; 

float G, AR, T,E,R, A,B,C, D, F,X, V,M1, Z,M2; 
scanf ( " % f % f 11 , &G, &AR) } 

printf ("The specific heat ratio = %f, The Area ratio = 
%f \n" , G, AR) ; 

switch ( i ) 

{ 

case 1: 

{ 

T=AR*AR; 

E=(G+1) / (G-l ) ; 

R=2/ (G+l) ; 

A=pow (R, E) ; 

B= ( l-T-2 *A) /A; 



C=B*B-4 ; 


X= (-B-sqrt (C) ) /2; 

Ml=sqrt (X) ; 

printf("The subsonic Mach number Ml = %f\n\n",Ml); 
break; 

} 

case 2 : 

{ 

T=AR*AR; 

E= (G+l ) / (G-l ) ; 

D=1 / (E-l ) ; 

Z=pow (T, D) ; 

V= (G-l ) / (G+l) ; 

R= ( G+ 1 ) / 2 ; 

A=pow (V, R) ; 

B= ( l-Z-2 *A) /A; 

C=B*B-4 ; 

X= (-B+sqrt (C) ) /2; 

M2=sqrt (X) ; 

printf("The super sonic Mach number M2 = %f\n\n",M2); 


break; 



Table 1: Data generated by using the above C-program 


Area ratio 

SUB SONIC MACH 

NUMBER 

SUPER SONIC MACH 

NUMBER 

1 

1 

1 

1.5 

0.424385 

1.788632 

2 

0.303366 

2 . 126523 

2 . 5 

0.238232 

2.371308 

3 

0.196688 

2.566362 

3.5 

0.167685 

2.726967 

4 

0.146225 

2.870822 

4.5 

0.129682 

2 . 995437 

5 

0.116523 

3 . 107255 
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